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^ I Abstract. Self-adjoint rank two commuting ordinary differential oper- 

C^l ■ ators are studied in this paper. Such operators with trigonometric, elliptic 

and rapid decay coefficients corresponding to hyperelliptic spectral curves 
are constructed. Some problems related to the Lame operator and rank two 
solutions of soliton equations are discussed. 
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1 Introduction 



^ I Almost all solutions of soliton equations obtained in the last decades with the help of 
I finite-gap theory are rank one solutions. This means that eigenfunctions of auxiliary 
linear operators form linear bundles over spectral curves. Meanwhile in the outstanding 
•/^ ■ papers of I.M. Krichever and S.P. Novikov [I], [2] rank I > 1 solutions are investigated 
CN . (i.e. the eigenfunctions form vector bundles of rank / over the spectral curves). One of 

the main difficulty to construct higher rank solutions (even numerically) is the problem 
of finding higher rank commuting operators. 

In [3] an example of rank two operators corresponding to nonsingular spectral curves 
■ of arbitrary genus is constructed. The operator 

l\ = {dl + a^x^ + a2x'^ + aix + ao)'^ + g{g + l)asx, 037^0 

commutes with a differential operator -^^4^+2 of order Ag + 2 [3J. The spectral curve 
is hyperelliptic curve of genus g. In this paper we further develop methods of [3] and 
construct self-adjoint operators of rank two with trigonometric, elliptic and rapid decay 
coefficients corresponding to hyperelliptic spectral curves. It is remarkable that there 
are higher rank operators with periodic and rapid decay coefficients expressed in terms 
of elementary functions (see Corollaries 1 and 2 below). The main results of this paper 
are Theorems 1 and 2. 



*This work was partially supported by grant 12-01-33058 from the Russian Foundation for Ba- 
sic Research; grant MD-5134. 2012.1 from the President of Russia; and a grant from Dmitri Zimin's 
"Dynasty" foundation. 
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Theorem 1 The operator 

l\ = {dl + a{P{x) + aof + aig2g{g + l)V{x), ^ 0, 
where V satisfies the equation 

{V'{x)f = g^Vix) + g{P{x) + ^o, 92 
commutes with an operator l\gj^2 of order Ag + 2. 
At go = 1, t/i = 0, g2 = -1 we have 
CoroUciry 1 The operator with periodic coefficients 

L\ — {dl + ai cos(x) + — aig{g + 1) cos(x), ai 
commutes with an operator of order Ag + 2. 

Let p{z) be the Weierstrass eUiptic function satisfying the equation 

{p'{z)f = V(^) + ^2P'(^) + gMz) + go. 

Theorem 2 The operator 

L\ = {dl + aip{x) + aof + Sip{x) + S2p'^{x), 

where 

ai = ^-2g'-2g, s, = ^g{g + l){16ao + 5g2), S2 = -ig{g + 2){g^ - 1), 

ao is an arbitrary constant, commutes with an operator of order Ag + 2. 
At g2 = 4a^ we have 

CoroUciry 2 The operator with rapid decay coefficients 

2 \ 2 2 4 

2 ctia \ Sia S2a 



where 



^ ^ cosh^(aa;) °/ cosh^(aa;) cosh'^(aa;)' ^ 



ai = ^-2g^- 2g, = ^g{g + l)(16ao + 200^), S2 = -Ag{g + 2)(/ - 1), 
CKo is an arbitrary constant, commutes with an operator of order Ag + 2. 

Equations of spectral curves for pairs L4, L\g_^_2 and L\, L\g_^2 s-^^ given in the Lemmas 
1 and 2. 



In the section 2 we recall the method of deformation of the Tyurin parameters. In 
the section 3 we prove Theorems 1 and 2. In the Appendix I the spectral curve and 
eigenf unctions of the Lame operator 

L2 = -dl + g{g + l)p{x) 

are found in some special form. In the Appendix II rank two solutions of soliton equa- 
tions are discussed. 

The author is grateful to A. P. Veselov for valuable discussions. 



2 Krichever— Novikov theory (method of deforma- 
tion of the Tyurin parameters) 

In this section we discuss some results on higher rank commuting ordinary differential 
operators and recall the method of deformation of Tyurin parameters (applications of 
this method for finding higher rank solutions of soliton equations see in PP,[2],[1])- 

Commutative rings of ordinary differential operators were classified by I.M. Krichever 
[5]. Common eigenvalues of commuting operators 

n—2 m—2 
i=0 i=0 

are parametriezed by points of the spectral curve F. If 

Ln^ = Z^, Lmlp = W^, (1) 

then {z, w) G F, where F is the smooth compactification of the curve defined by the equa- 
tion R{z, w) = 0, R is the Burchnall-Chaundy [6J polynomial such that R{Ln, L^) = 0. 
The rank of the pair {Ln, Lm) is called the dimension of the space of common eigenfunc- 
tions ([1]) for fixed P = {z,w) & T in general position. Eigenfunctions (Baker-Akhiezer 
function) and coefficients of rank one operators can be expressed via theta-function of 
F [Zj. Operators of rank two corresponding to the elliptic spectral curves were found 
by I.M. Krichever and S.P. Novikov yy, operators of rank three with the same spectral 
curve were found by O.I. Mokhov [1]. These operators and some related problems were 
studied in the papers [10]-[8]. Some examples of operators of rank 2 and 3 with spectral 
curves of genus 2,3,4 were found in [T8] [21]. 

The case of rank 1 the Baker-Akhiezer function ip{x, P) has zeros divisor 7(x) = 
7i(x) -|- ■ ■ ■ + 7g(a;) where g is genus of F. The x-trajectory of 7 in the Jacobi variety of 
F is a straight line. In the case of higher rank corresponding x-dynamics in the moduli 
space of vector bundles over F is very complicated. Let me recall spectral properties 
of ip = {ipi, . . . at Z > 1 [5]. The Baker-Akhiezer function ip has Ig simple poles 
Pi, . . . ,Pig with the properties 

Resp^ipj = VijReSp^ipi, I <i < Ig, 1 < j < / - 1. 
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The set 



{iPj,Vik)} 



(2) 



is called Tyurin parameters, which define a holomorphic bundle of rank / over F. There 
is a point g G F where ip has essential singularity of the form 



where eo = (l,0,...,0),6(a;) = (CHa;) 
^ satisfies the equation 



ii\{x))i is a local parameter, the matrix 



,s=0 



cia; 
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\k + Uo{x) uji{x) ui2{x) 








\ 
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(^i-2ix) / 



The function can not be found explicitly. To find operators one can use the 
following method of deformation of Tyurin parameters. Let us consider for simplicity 
the case of rank two. Then for P G F the space of common eigenf unctions satisfies the 
second order differential equation 

tlj"-Xi{x,P)tlj'-Xo{x,P)tp = 0. 
The operator — xA — Xo is a common right divisor of L„ — z and — w. Functions 



Xo, Xi are rational functions on F, with poles at pi{x) 
pole at q, Xi has zero at q 



,P2gix) and Xo has additional 



Xo = k + O{l), xi = Oil/k). 
Let k — 'ji{x) be a local parameter near Pi{x), then 

Xo{x, P) = -y^'^^f'^^ + + 0{k - 7.(a:)), 

~ liyx) 

Xiix, P) = . + di,^ix) + Oik - 7.(x)), 

where di^, di^i, ■ji{x) , Vi^o{x) satisfy the Krichever equations [5] 

di,o{x) = ViQ^x) + Vi^o{x)di^i{x) - v[q{x), i = l,...,2g. 



(3) 



The last equations define deformations of Tyurin parameters and {pi{x) , Vifl{x)) coin- 
cides with ([2]) at some x = Xq. I.M. Krichever and S.P. Novikov found Xo^Xi and 
corresponding operators at g = 1,1 = 2. 
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Theorem 3 (I.M. Krichever, S.P. Novikov JT]) The operator 

Lkn = {dl + u)^ + 2c^(p(72) - p(7i))5x + (cx(p(72) - p{li)))x - ^(72) - p{li), 

liix) =70 + c(x), 72 (x) = 7o - c(x), 
u{x) = -i^ + 1% + 2$(7i,72)c. - + 4($,(7o + c,7o - c) - $^(71,72)), 

$(71, 72) = C(72 - 7i) + C(7i) - ((72), 

where c{x) is an arbitrary smooth function, 70 is a constant, commutes with an operator 
of order six. 

The very interesting examples of the Krichever-Novikov operators are Dixmier op- 
erators. 

Theorem 4 (J. Dixmier 12^ ) The operators 

\dx?^ ) 

commute, herewith 



The action of the automorphisms group of the first Weyl algebra on the Dixmier op- 
erators and on the operators L\,Llg_^2 were studied in [23]. The conditions when Lkn 
has rational coefficients were found in [H]. 

Theorem 5 (P.G. Grinevich 19]) The operator Lkn has rational coefficients if and 
only if 

°° dt 

where q{x) is a rational function. // 70 = and q{x) = x, then the operator Lkn 
coincides with the Dixmier operator L^. 

The conditions when Lkn is self-adjoint were found in p^ . 

Theorem 6 (P.G. Grinevich, S.P. Novikov fTUf ) Operator Lkn is formally self-adjoint 
if and only if ^(71) = ^(72). 

In [3] the operators of rank two corresponding to the hyperelliptic spectral curve 

T:w^ = Fg{z) = z"^^^ + c^gz"' + ■ ■ ■ + Cq 




5 



were considered 

L4%i) = zip, = wip, P = {z, w) e r. 



Theorem 7 (M. JE]) The operator L4 is self-adjoint if and only if 

Xi{x,P) = Xi{x,cr{P)), (4) 

where a is the involution 

a{z, w) = {z, —w). 

Let us assume that xi is invariant under a, then L4 can be represented as 

L4 = {dl + V{x)f + W{x) 

and xo, Xi have the form 

Theorem 8 (M. J^]) Functions XoiXi have the form 

IQ" w _Q' 

where Q is a polynomial in z of degree g with coefficients depending on x 

Q{x, z) = {z- 7i(x)) ...{z- 79(3;)). (5) 

The polynomial Q satisfies the equation 

AFg = A{z - W)Q^ - AViQ'y + (Q")' - 2Q'Q^''^ + 2Q{2V'Q' + AVQ" + g(^)), (6) 

where Q' , Q" , Q^'^^ mean d^Q, d^Q, d^Q. 

Corollary 3 The function Q satisfies the linear equation 

C,Q = {d'^ + 2{V, dl) +2{z-W- V", d,)) g = 0, (7) 

where {A, B) = AB + BA, the operator £5 is skew- symmetric. Potentials V, W have 
the form 

\Q"f -200^^ -AFg 
4(Q0 

W = -2(71 + . . . + 7j - C2g. 

The functions 7i(x), . . . , 7g(a;) satisfy the equations 



iQ'r~2Q'Q(^^-AFgiz) _ {Q"r-2Q'Q(-'^-AFg{z) 

4(Q')' 4(Q')' ■ 
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So, Theorem 8 and Corollary 3 show that the spectral theory of the operator L4 
is similar to the spectral theory of the finite-gap one dimensional Schrodinger operator 
(see 121]) 

{—dl + u{x))ip = zip. 

To the Schrodinger operator also corresponds a skew-symmetric operator but of order 
three [21] (see also [2S]) 

CsQ= {dl + 2{z-u,d,))Q = 0, 
and Q satisfies the nonlinear equation 

4Fg{z) = 4{z - u)Q^ - {Q'y + QQ", (10) 
The potential u has the form 

U = -2(71 H h 7g) - C2g. 

Equations ([9]) are analogues of the Dubrovin equations 

Using the equation V.N. Davletshina [26] proved that the operator 

{dl + aie"" + aoY + g{g + l)aie'', cti 7^ 
commutes with an operator of order Ag + 2. 

3 Proof of Theorems 
3.1 Theorem 2 

We construct the polynomial satisfying the equation ([71) for 

= aiV{x) + ao, = aig2g{g + l)V{x). 

Let 

Q\x,z) = Al{z)V^ix) + --- + Al{z). 
Let us substitute Q'' into ([7]) and use the identities 

V"{x) = g^Vix) + |, V"'{x) = g^V'ix), 
V^'\x) = glV{x) + P(^)(x) = ^^^(x). 
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We get 

1 



4 

where 

si si si 

+8gog2{s{s + 4) + 5) + 2^+2^^^^(2s + 3){Aa,go + g,{Aao + g2{2s{s + 3) + 5))) + 

4As+i{s + + l)\4aigi + gii^a^gi + ^2(4ao + g2{s + 1)' + 4z) + 
SAs{2s + l)aig2{s{s + 1) - g{g + 1)) = 0, 
we assume that = at s < and s > g. Hence from Ps = 0,0 < s < g we have 

As = ^TTT—T^ — rr—T\ — r—u\ ( 4^+5^o^^^^ + 4^+4(70^71^^^^ (2s + 5) 

8{2s + l)aig2{s{s + 1) - g{g + 1)) \ si si 

+As+3^^^^{lQao9o + ^7i(2s + 3){2s + 5) + 8(70^72(^(5 + 4) + 5)) + 
si 

+AAs+i{s + l)((s + l)^(4ai^i + ^2(4ai^i + ^2(4ao + g2{s + if + Az)) . (11) 

Let 

, _ (1 ■ 3 (2(7 + 1K(2- (7(^ + 1))... ((^-1^-^(^7 + 1)) .... 
- ^\ • 

Then Q'' has the form (jS]) and satisfies ([7]) and ([6]). Theorem 2 is proved. 
Let us find the spectral curve. 

Lemma 1 The spectral curve of L\,L\g_^_2 is given by the equation 

= ^{AAlz+Ao{16aogoA2+A8glAi+3QgogiAs+3gfA2+l6gog2A2+Ai{{25~8g{g+l))go 

+^i(4ao + ^2))) - 4ao^o^? + ^(4^o^ + 9iA? - 2^o^i(6^o^ + S^i^ + ^2^1)), 
where Aj are defined in 07]) ./ fl^) . 

Proof. Let xq be a zero of V{x). Since in the right hand side of ([6]) there are only 
derivatives of order not grater than four, then we have 

4F, = ^{4{z - W)Q' - 4V{Q'f + {Q")' - 2Q'Q^'^ + 2Q{2V'Q' + 4VQ" + Q'^'^))U=.„ 
where 

Q^ = A\iz)V\x) + --- + Aliz). 

Using 

iV\x)f = go, V"{x) = ^, V^'\x) = g2V'{x), V^'\x) = ^-^ 
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we get the equation. Lemma 1 is proved. 
Let 

H'^ = d1 + «! cos(x) + tto, L\ = (H^Y — aig{g + 1) cos(x) 
and we shall use the notation 

<A,B>= AB + BA. 

Examples: 

1) ^7 = 1 

L\ = {H^f + - < 1 - 4ao - 12ai cos(x), > +ai cos{x), 
8 

Ff{z) = z^ + {l- 2ao)z^ + -^(1 - 8ao + 16a^ - I6al)z + ^. 

2 lb 4 

2) g = 2, let ao = 

L\q = [H^f + 1 < II _ I5ai cos(x), {H^f > -—ai < cos(x), {H^Y > + 

1 3 

- < 1 + 27^1 - 60ai cos(x) + Abaj cos(2a;), H'^ > +-ai(5 cos(a;) - 12ai), 

2 2 

F^{z) = z^ + —z^ + —(321 - 336al)z^ + -{34 - 531al)z^+ 
2 16 4 

(1 - 189a? + 108at)z + 24a? + 513at 

3.2 Theorem 3 

Let 

Q\x,z) = Aliz)p\x) + --- + Aliz). 
Let us substitute into ([7]). Using the identities 

p"{x) = (6p2(x) + g2p{x) + y), <9^p(x) = {12p{x) + g2)p'{x), 

dtp{x) = {120p%x) + 30g2p\x) + {18g, + gl)p{x) + ^ + 12go), 
dlpix) = {360p\x) + 60g2pix) + 18g, + g!)p'ix) 

we get 

C,Q = p'{x)i(3g+,p^^\x) + ■■■ + p{x)) = 0, 
where (3s+i = is equivalent to 

16A,{s + 1)((6 + 8ai)s + (19 + 4:ai)s^ + I6s^ + 4/ - ^2) + 

8A,+i(2s + 3)(4ao(s + l)(s + 2) + ^2(5 + l)(s + 2){2s'^ + 6s + ai + 5) - + 
4A,+2(4ao^2(s + 3)^ + gl{s + 2)^ + 2g,{s + 2f{A{s + 2f + 2ai + 5) + 4(s + 2)z) + 
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2A,+3(s + 3)(s + 2)(2s + 5)(4^o(«i + 2(5^ + 5s + 9)) + ^i(4ao + ^2(25^ + 10s + 13))) + 
^+4(5 + 2)(s + 3)(s + 4)(16ao^o + 8^0^2(5^ + 6s + 10) + gl{As^ + 24s + 35)) + 

AAs+^gogi{s + 2)(s + 3)(s + 4)(s + 5)(2s + 7) + 4^+6^7o|j^ = 0- 

At s = (? and s = g — 1 this formula has the forms 

Ag{g + 1)((6 + 8ai)(7 + (19 + Aa^)g^ + IQg^ + Ag"" - s^) = 0, (13) 

Ag{2g + l){AaQg{g + 1) + (?((? + 1)(1 + ai + 2(7 + 2g^)g2 - si) + 

2A,_i(7(l - hg^ + 4(7^ + 4ai((7^ - 1) - ^2) = 0. (14) 

At 

Si = ^g{g + l){lQa^ + hg2), S2 = -Ag{g + 2){g^ - I), a^ = ^-2g-2g^ 



equations (jT3l) . ([T4|) are vanished. From = we get 

A,_i(4ao(8(? - 4) + {IQg^ - 2Ag^ + l^g - 5)^72) 



64(2(?2_ 3(7 + 1) 

^g(8^'^?i - g^gl -g\llgi + Aa^g^) - Az) 
64(2(72-3(7+1) 

Thus from the equations /3j = 0, z > we express v4j_i through Ag^i and Ag. The last 
equation /3o = takes the form 

^gP0(^)+ Vl^l(^) =0, 

where Po{z) and -Pi (-2) are some polynomials. Let 

Ag = -Pl{z), Vl=P0(^). 

So, we find polynomial in z solution of the equation (1). Theorem 3 is proved. 

The equation w'^ = Fg{z) of the spectral curve of L\, L\g_^_2 is given in the following 
lemma. 

Lemma 2 The spectral curve of L\^L\gj^2 given by the equation 



w 



2 - - [AAlz + Ao(4«i Ai^o + 48^2 + 36^3^0^! + 'iA^gl + Aa^iAA^g^ + A^g^) + 16A2^o^2+ 
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+Ai^i^2) - Aa^Alg^ + i(4A2^o + A^g^f - 2Ai^o(6A3^o + 3^2^! + A^g^)), 
where Aj are defined above. 

The proof of the Lemma 2 is the same as the proof of the Lemma 1. 
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Let 

= dl + aip{x) + ao- 

Examples: 

For the simplification of the formulas we put g2 = ^- 

3) (7 = 1 

4 = {H'f + l<gi + 2aop{x), > +2aogi + 32aop\x), 

o 

Ft = z'+ ^-fz' + {Qaogo + Hdi + ^5?)^ + Hdl + ^^^^ - l^o^ho. 

4) (7 = 2,^1 = 

4p = {H^f+ < 30aop(a;)-12p2(x), (if^)^ > +15 < 16^o-12aop^(x)+16p^(x), (i/^)^ > 

+9 < aogo - 1850^op(a;) + 80alp'^{x) - 480aop^ - 11520p^(x), > - 
36{28algo - llSlao^oP - Q9A6gop^{x) + 160a^p^(x) - 7520aop\x) - 30208p^(a;)), 
F^{z) =z^ - 387aogoZ^ + 27^o(16aj^ + 139^0)^^ + 12636a^c/^z - 243aog^{64al + 637^o)- 

4 Appendix I. Spectral curve and eigenfunctions of 
the Lame operator 

S.P. Novikov proved [27] that if the periodic Schrodinger operator 

-dl + u{x) 

commute with an operator of odd order then it has finite number of gaps in its spectrum 
(the inverse statement were proved by B.A. Dubrovin [28]). Such operator is called 
finite-gap operator. Eigenfunctions of the finite-gap operator are expressed in terms 
of theta-function of the hyperelliptic spectral curve via Its formula [29]. The famous 
important finite-gap operator is the Lame operator 

L2 = -dl+gig + l)pix). 

It was introduced in 1837. Hermite and Halphen used the following ansatz for eigen- 
functions of the Lame operator 

= 

<;(9+i)/2 

^= aj{z,aj)^{x,aj), 

where ^ ^ 

a[a)a[x) 
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aj{z,aj) is some function. With the help of this ansatz Hermite and Halphen studied 
the cases g = 2,3. I.M. Krichever [30] used some modification of the ansatz to find 
elhptic solutions of the Kadomtsev-Petviashvili equation (see also [31]). M.-P. Grosset 
and A. P. Veselov [32] proved that coefficients of the equations for spectral curve of the 
Lame operator can be found recurrently with the help of elliptic Bernoulli polynomials. 

In the next theorem we find the equation of the spectral curve in the explicit form. 
For simplicity of the formulas we assume in this section g2 = 0, i.e. the Weierstrass 
function satisfies the equation 

{p'{z)f = 4:p%z)+gip{z)+go. 

Let us introduce the following polynomials in z with the help of recurrent formulas. Let 

As = at s > g and s < 0, 

. _ i<^{g^ + g)..{g^ + 9-{s + l)s)..{g^ + 9 - {g - l)g) 

^ _ (g + l){SAs+iz + As+2gi{s + 2){2s + 3) + 2As+^go{s + 2)(g + 3)) 

A{2s + l){g^ + g-s{s + l)) 

We also introduce the following function 

Q = A,p3{x) + --- + A^. (15) 



Theorem 9 The spectral curve of the Lame operator is given by the equation 

= ^{4Alz + Ao{4A2go + A^g,) - Ajgo). 



The eigenfunctions of the Lame operator are 
Proof. Our observation is that the solution of the equation (ITOl) for the Lame operator 



is (ITSll . The proof is the same as proof of Theorems 1 and 2. To find eigenfunctions we 
observe that the eigenfunctions satisfies the equation 

Qx w 
2%Q Q 

It follows from 

-dl + g{g + l)p(x) = (-9^ - ixo){d - ixo)- 

Theorem 9 is proved. 

The last theorem gives very effective way to find spectral curve of the Lame operator. 
Examples: 
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5) g = l 

6) g = 2 

7) (7 = 3 



r 3 I fi'i 9o 
F, = Z "r-Z--. 



5 21^1 3 27 2 27 2 81 
^2 = + + ^g^z" + —g\z + —(70^71. 



7 63 , 297 4 4185 o ^ 18225 o 3375 o 

^3 = + ^9x'^' + ^^70^^ + ^TT^i^ + ^^^?o^7i^' + ^7r(27^?o + 9\)- 
/ z 16 8 16 

5 Appendix II. Rank two self-adjoint operators and 
evolution soliton equations 

Let us consider the following system of evolution equations 

Vt = ^{QVV, + 6W, + V,,,), Wt = ^{-3VW, - W,,,). (16) 
This system is equivalent to the commutativity condition 

where 

L, = {dl + V{x,t)f + W{x,t), ^3 = 5' + \v{x,t)d, + -V^{x,t). 

Following to the papers of I.M. Krichever and S.P. Novikov [H |2] we call the solution 
V{x, t), W{t, x) the solution of rank I if for every t the operator L4 is included in a pair 
of commuting operators of rank / (hire Z = 1 or Z = 2). Rank 1 solutions of (|T6l) were 
found by V.G. Drinfeld and V.V. Sokolov [33j. It is a very natural question what is the 
evolution equation for Q under (1161) ? 

Theorem 10 (V.N. Davletshina, M.) The polynomial 

Q = (z- 7i(a;, t))...{z- -fg{x, t)) 

satisfies the equation 

Qt = ^i-SVQ,-Q,,,). (17) 



We shall give the proof of the Theorems 10 and 11 and related statements in an 
another publication. 

The equation ([T7|) gives a symmetry of the equation (Q. At g = 1 the formula ([8]) 
gives 

^ -16Fi(|(-C2 - W)) + W^, - 2W,W,,, 
13 



where 

nP' = Fi{z) = + C2Z^ + ciz + Co, 



is the equation of the spectral curve, and the equation flTTI) is reduced to the famous 
Krichever-Novikov (KN) equation 

_ 48Fi(i(-C2 - W)) - W^, + 2WJ¥,,^ 
Wt — 



8W^ 

KN equation plays an important role in the theory of rank two solutions of the Kadomtsev- 
Petviashvili (KP) equation at g = 1. It would be very interesting to study the equation 
( !T7|) in the framework of rank two solutions of KP ai g > 1. 



Theorem 11 (V.N. Davletshina, M.) The system [T^) has the following rank two solu- 
tion corresponding to the elliptic spectral curve 

20 

W{x, t) = -AOp{bt + x)^ - —(86 + 7g2)pibt + x) + c, 



Vix,t) = -10pibt + x)-^-^, 



where b,c — const. 



Let us numerically study the Cauchy problem for f|T6l) with initial data from Corol- 
lary 2 and with the rapid decay boundary conditions 

V{x, 0) = — - + ao, 

cosh [ax) 

2 4 

W{x, 0) = 



cosh^(ax) cosh^(ax)' 

\4(±oo,t) = V;,(±oo,t) = W,i±oo,t) = W,,{±oo,t) = 0. 

The behavior of the solution looks like (yf-soliton equation of the Korteveg-de Vris equa- 
tion (see Fig. 1 — 12). 

It is an interesting problem to find exact solution of the Cauchy problem. It gives 
soliton deformations of the rank two self-adjoint commuting differential operators. 



Figure 1: one-soliton solution, g = l,a = l/2,aQ = 0,t = 

14 



w 

0.2 - 




-O.S - 

Figure 2: one-soliton solution, (7 = 1, a = 1/2, ao = 0, t = 



-40 -20 



20 40 



Figure 3: two-soliton solution, g = 2,a = 1/2, = 0, t = —15 



-40 -20 



20 40 



Figure 4: two-soliton solution, g = 2,a = l/2,ao = 0,t = —3 



Figure 5: two-soliton solution, g = 2,a = l/2,ao = 0,t = —1.5 



Figure 6: two-soliton solution, g = 2,a = 1/2, = 0, t = 1 
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Figure 7: two-soliton solution, g = 2,a = l/2,ao = 0,t = 5 




Figure 8: three-soliton solution, g = 3,a = 1/4, ao = 0, t = —30 



Figure 9: three-soliton solution, g = 3,a = 1/4, ao = 0, t = —15 



Figure 10: three-soliton solution, (7 = 3, a = 1/4, ao = 0, t = —7 



Figure 11: three-soliton solution, g = 3,a = 1/4, ao = 0,t = 4 
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Figure 12: three-soliton solution, g = 3,a = 1/4, ao = 0,t = 20 
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